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ABSTRACT 

We investigate the star formation and chemical evolution in the early universe by considering the 
merging history of the Galaxy in the A cold dark matter scenario according to the extended Press- 
Schechter theory. We give some possible constraints from comparisons with observation of extremely 
metal-poor (EMP) stars, made available by the recent large-scaled surveys and by the follow-up high- 
resolution spectroscopy. We demonstrate that (1) The hierarchical structure formation can explain 
the characteristics of the observed metallicity distribution function (MDF) including a break around 
[Fe/H] = —4. (2) A high mass initial mass function (IMF) of peak mass ~ 10 Mq with the contribution 
of binaries, derived from the statistics of carbon enhanced EMP stars (Komiya et al. 2007, ApJ, 658, 
367), predicts the frequency of low-mass survivors consistent with the number of EMP stars observed 
for —4 < [Fe/H] < —2.5. (3) The stars formed from primordial gas before the first supernova 
(SN) explosions in their host mini-halos are assigned to the hyper metal-poor (HMP) stars with 
[Fe/H] ~ —5. (4) There is no indication of significant changes in the IMF and the binary contribution 
at metallicity —4 > [Fe/H] > —2.5, or even larger as long as the field stars of Galactic halo are 
concerned.. 

We further study the effects of the surface pollution through the accretion of interstellar matter 
(ISM) along the chemical and dynamical evolution of the Galaxy for low-mass population HI and 
EMP survivors. Because of shallower potential of smaller halos, the accretion of ISM in the mini- 
halos in which these stars were born dominates the surface metal pollution. This can account for the 
surface iron abundances as observed for the HMP stars if the cooling and concentration of gas in their 
birth mini-halos is taken into account. 

We also study the feedback effect from the very massive population HI stars. The metal pre- 
pollution by pair-instability SNe (PISNe) is shown to be compatible with the observed lack of their 
nucleosynthetic signatures when some positive feedback on gas cooling works and changes IMF from 
being very massive to being high mass. 
Subject headings: Galaxy: formation - Galaxy: evolution - stars: abundances - stars: Population II 



1. INTRODUCTION 

In the early universe, the first stars should be formed 
from the primordial matter without metal elements (e.g. 
Barkana & Loeb 2001; Bromm & Larson 2004, and ref- 
erences therein). These stars formed without metals are 
referred to as Population HI stars. Despite longstand- 
ing efforts, however, direct proof of their existence has 
been elusive. The first stars illuminate the universe and 
enrich surrounding matter with metals and the subse- 
quent generations of stars with low metallicity should 
be formed under their influence. At the early stages of 
formation and evolution of galaxies, the feedback from 
stars is thought to play an important role. Since the era 
of the first stars and galaxy formation has occurred in 
the Dark Age of the universe, however, it is difficult to 
investigate them by direct observation with current in- 
struments. One promising observational probe into this 
era is the extremely metal-poor (EMP) stars, which were 
stars formed in the early stages of the Galactic chemical 
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evolution. The observations of EMP stars are expected 
to be connected with the first stars and Galaxy forma- 
tion. 

Over the past decade, thousands of very metal-poor 
stars below [Fe/H] ~ —2 have been identified in the 
Galactic halo by two large-scaled surveys (HK sur- 
vey. Beers et al. 1992: Hamburg/ESO [HES] survey, 
Christlieb et al. 2001). The follow-up observations with 
high-dispersion spectroscopy (HDS) using large tele- 
scopes yield details of hundreds of these star's surface 
abundances (e.g., see the review by Beers & Christlieb 
2005). Thanks to these observations, now we are able 
to discuss statistical features of these stars (see, e.g., 
SAGA database; Suda et al. 2008, 2009). It is partic- 
ularly worth noting that the most iron-deficient objects 
known to date are identified among the HES sample with 
much smaller metallicity than other EMP stars, i.e., with 
[Fe/H] = -5.6, -5.4 and -4.8 (HE1327-2326, Frebel et 
al. 2005: HE0107-5240, Christlieb et al. 2002: HE0557- 
4840, Norris et al. 2007). We refer to these three stars as 
hyper metal-poor (HMP) stars, distinct from EMP stars, 
defined as stars with [Fe/H] < —2.5. 

The EMP stars display some features distinct from the 
Population II stars and Pop. HI stars, as summarized in 
Section 2.1. In previous studies (Komiya et al. 2007, 
2009a, , hereafter Paper I and Paper II, respectively) we 
propose that they compose a peculiar stellar population 
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formed with a high-mass initial m.ass function (IMF). Wc 
also argue that most of observed EMP stars were formed 
as the secondary members in binaries. 

From the observations of EMP stars in the Galac- 
tic halo, the scarcity of stars below [Fe/H] ~ —4 has 
emerged and has been investigated from the perspective 
of chemical evolution (Prantzos 2003; Karlsson 2005). 
Among - 150 stars with [Fe/H] < -3 with HDS ob- 
servations, only three HMP stars are discovered be- 
low [Fe/H] ^ —4.8 and no stars are identified between 
—4.8 < [Fe/H] < —4, although a few stars are found 
slightly below the [Fe/H] — —4. In Paper II, we pro- 
pose a new scenario for this scarcity that the stars are 
divided into two groups according to whether they were 
formed before and after the parent mini-halos were pol- 
luted by their own first supernovae (SNe). HMP stars 
belong to the former group while the stars of the latter 
group constitute EMP stars. We summarize this scenario 
in Section 2.2. 

There have been investigations on the origins of the 
metal elements in the HMP stars. All of three HMP 
stars, known to date, show large carbon enhancement. 
Suda et al. (2004) propose a binary scenario to assert 
that HMP stars can be the survivors of the low-mass 
Pop. HI stars without pristine metal elements. In this 
scenario, their surface metal elements originate from af- 
terbirth pollution by the accretion of interstellar mat- 
ter (ISM) and by the wind accretion of envelope matter 
ejected by the binary companions. Umeda & Nomoto 
(2003) and Iwamoto et al. (2005) propose that HMP 
stars are formed of the gas, enriched with metals by the 
ejecta of peculiar faint Pop. HI SN. Limongi et al. (2003) 
propose that a combination with two Pop. HI SNe can 
reproduce the abundance pattern peculiar to HMP stars. 
Maynet et al. (2006) suggest that light metal elements on 
HMP stars are brought by the stellar wind from rapidly 
rotating Pop. HI massive stars. For the binary scenario, 
Nishimura et al. (2009) investigate the AGB nucleosyn- 
thesis under the dearth of pristine metals to show that 
the abundance patterns of light metal elements from C 
through Al can be reproduced in terms of the accretion 
of envelope matter from AGB companion in binaries. In 
this case, the origin of the iron group elements remains to 
be discussed because binary mass transfer or stellar wind 
cannot provide these elements. As for the SN scenario, 
there remains a problem of forming the next generation 
of stars with an appropriate mixture of SN ejecta and 
ISM, and in particular, for faint SNe with the peculiar 
abundances, the stars have to be formed with the effi- 
cient use of a relatively small amount of SN ejecta. It is 
necessary, therefore, to discuss the origin of iron on the 
surface of HMP stars in the proper context of structure 
formation of the Galaxy. 

From the theoretical viewpoint, many calculations of 
the formation of first stars in the early universe have 
been carried out (e.g. Bromm et al. 1999; Abel et al. 
2002; Omukai & Para 2003; Tan & McKee 2004; O'Shea 
& Norman 2007; Yoshida et al. 2008b). These studies 
suggest that one very massive star with M > 100 M© 
may be formed in the primordial mini-halos with mass 
10^ Mq, although typical masses or the mass range of 
the first stars are yet to be properly decided, especially 
in relation to the role of angular momentum in the frag- 
mentation (Clark et al. 2008; Machida et al. 2008). On 



the other hand, the evolution of Pop. HI and EMP stars 
have been investigated rather well, and in particular, it 
is shown that very massive stars with M ~ 200 Mq ex- 
plode as pair-instability SNe (PISNe) with huge explo- 
sion energy and a large amount of metal yield (Umeda 
& Nomoto 2002; Heger & Woosley 2002). The effects of 
these energetic SNe on the subsequent star formation and 
galactic chemical evolution are not yet understood well, 
however (e.g. Machida et al. 2005; Kitayama & Yoshida 
2005; Greif et al. 2008). The conditions for the forma- 
tion of the first low-mass star are subject to controversy 
as well. Some studies propose that metal enrichment by 
massive stars accelerates the cooling of molecular clouds 
to lower the Jeans mass and to enable the low-mass star 
formation (Bromm & Loeb 2003; Omukai et al. 2005). 
Other studies propose that ionization photons emitted 
from first stars accelerate the formation of H2 molecules 
as main coolant in the primordial gas (Ricotti et al. 2002; 
Yoshida et al. 2008a), which can trigger the first forma- 
tion of low- mass stars. There remains the possibility that 
low-mass stars can be formed as the first objects (Naka- 
mura & Umemura 2001; Clark et al. 2008). 

In this paper, we focus on the origin and nature of 
HMP and EMP stars and discuss the first stars, the early 
chemical evolution, and the star formation history of the 
Galaxy in the context of the hierarchical structure for- 
mation scenario. We have built the merger trees of the 
Galaxy using a semi-analytic method and follow the his- 
tory of star formation and chemical enrichment along the 
trees (Komiya et al. 2009b, hereafter, Paper HI). We reg- 
ister all the individual Pop. HI and EMP stars and trace 
their formation and evolution. 

Firstly, wc derive the MDFs and compare them with 
those of the observational counterparts in the Galactic 
halo. We adopt the high mass IMF, based on the com- 
parisons of stellar models with observations, and take 
account of the binary contribution according to Paper I 
and Paper II. We also illustrate the chemical enrichment 
history in the hierarchical structure formation scenario 
and discuss the dependence of MDF on the assumptions 
about feedback from first stars and star formation rate 
in the early universe. In our study, we take account not 
only of the shape of the MDF but also of the total num- 
ber of EMP and HMP stars. 

Secondly, on the basis of this scenario, we study the 
effect of surface metal enrichment by accreting the ISM 
and attempt to predict the observational abundances of 
the first low-mass stars. Unlike the Pop. I and Pop. II 
stars, the ISM accretion can affect the surface metal 
abundance of Pop. HI and EMP stars because of their 
very low pristine metal abundances (Yoshii 1981; Iben 
1983). Furthermore, in small halos, in which Pop HI and 
EMP stars are born, the accretion rates of ISM gas are 
thought to have been much higher as compared with the 
those in the present Galactic halo because of much shal- 
lower gravitational potential. In order to estimate the 
ISM accretion rate and the change in the surface abun- 
dances of halo stars in the early universe, therefore, we 
should follow the merging history of halo as well as the 
star formation and chemical evolution. There are some 
previous studies of the ISM accretion (Iben 1983; Frebel 
et al. 2009) onto Pop. HI stars but none of them take 
account of the dynamical and chemical evolution self- 
consistently. In this paper, we present the calculations 
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of ISM accretion taking into account the merging history 
and the metal enrichment history of the Galaxy to ex- 
plore the dependence of accretion rates on the dynamics 
of gas and stars. We demonstrate that the ISM accretion 
provides a reasonable interpretation of the observed iron 
abundance of HMP stars. 

Thirdly we study possible metal pre-poUution by ener- 
getic SNe and draw observational constraints on the first 
stars and feedback from them. If PISNe or hypernovae 
(HNe) with the explosion energy and iron yield much 
larger than those of normal core collapse SNe (Heger & 
Woosley 2002; Umeda & Nomoto 2002) occur in the low- 
mass halos, the blast wave may sweep out the gas from 
their host halos (Bromm et al. 2003; Kitayama & Yoshida 
2005; Machida et al. 2005). As a result, their metal ejecta 
is blown off from the host halos to pollute the inter-halo 
gas. We calculate the MDFs with such effect due to 
PISNe according to the hierarchical chemical evolution 
program and compare them with the observed HMP and 
EMP stars. 

A semi-analytic approach to hierarchical structure for- 
mation has been adopted by Tumlinson (2006) and by 
Salvadori et al. (2006), who have constructed the merger 
trees of the Galaxy and investigated inhomogeneous 
Galactic chemical evolution models. Tumlinson (2007a) 
considered an IMF dependent on the temperature of the 
cosmic microwave background, and discussed the carbon- 
rich EMP stars also by taking into account the contribu- 
tion of binaries. However, he did not relate hierarchical 
nature to the scarcity of the HMP stars. We address the 
difference between the HMP and EMP stars in the con- 
text of hierarchical galaxy formation. Additionally, we 
discuss the total number of EMP stars as well as shape 
of the MDF. Furthermore, none of these studies is ca- 
pable of explaining the origin of the extraordinarily low 
iron abundance of HMP stars. We examine two possible 
origins, i.e., the surface pollution by ISM accretion and 
the pre-pollution of intra-Galactic matter (IGM) by halo 
blown off, for the small amount of iron in HMP stars. 

The paper is organized as follows. We review the ob- 
servations and elaborate our scenario in the next section 
and describe the computational method of star forma- 
tion and evolution of Galaxy, in Section 3. The following 
three sections are devoted to presenting the results of our 
computations for the halo formation models, the surface 
pollution due to the ISM accretion, and the effects of pre- 
pollution of mini-halos by energetic PISNe, respectively. 
We discuss the formation of the first stars and HMP stars 
in Section 7 and a summary of the conclusion follows in 
Section 8. 

2. OBSERVATION AND SCENARIO PRELIMINARIES 

In this paper, we use the following terminology for the 
classification of metal-poor stars. 

EMP stars: Stars with metallicity — oo < [Fe/H] < 
—2.5 different from those with of — 4[Fc/H] < —3 in the 
literature like Beers & Christlieb (2005) because the 
stars with [Fe/H] < — 2.5 show some theoretical and ob- 
servational peculiarities, distinct from the stars of larger 
metallicity, as described in Section 2.1. Stars formed 
with Z = are classified into the other population. 
Pop. HI stars. 



EMP population: Population of stars with metallicity 
— oo < [Fe/H] < — 2.5. It is the mother population of 
low-mass EMP survivors, and involves the stars with 
mass larger than ~ 0.8 M© that have ended the nuclear 
burning stages. 

EMP survivors: EMP stars with nuclear burning 
still going on. They are the low-mass members of the 
EMP population with mass M < 0.8 Mq. Most of these 
should be the secondary companions of binaries under 
our assumed high- mass IMF (see §2.1 and Paper I). 
Among the EMP stars born in the early universe, it is 
only EMP survivors that are currently observed in the 
Galactic halo. 



HMP stars: 
[Fe/H] < -4.5. 
to this category. 



Stars observed with mctallcity 
Currently, three stars are assigned 



Pop. Ill stars: Stars formed in primordial gas totally 
devoid of metals in their interior. In addition to the first 
stars, some may be formed under the radiative feedback 
effect from the first massive stars. Low-mass Pop. HI 
stars are, if formed, still shining as with the nuclear 
burning today, and are referred as "Pop. Ill survivors ". 
Pop. HI stars may suffer from the surface pollution by 
ISM accretion after their birth and appear enriched with 
finite metal abundances on their surface today. We refer 
to these stars, which are formed as Pop. HI stars but no 
longer Z = hy surface pollution, as "polluted Pop. Ill 
stars". In section 5, we discuss the possibility that HMP 
stars are polluted Pop. HI stars. 

2.1. Peculiarities of Stars with [Fc/H] < —2.5 

EMP stars are formed of gas that contains small but 
finite abundances of metal elements, synthesized in stars 
of previous generations, and are distinct from Pop. HI 
stars. These stars are also distinguished from Pop. II 
stars by peculiar features, as summarized below. 

Theoretically, the evolution of low- and intermediate- 
mass stars of metallicity [Fe/H] < —2.5 differs from that 
of Pop. II stars (Fujimoto et al. 1990, 2000). A mixing 
mechanism called He-flash driven deep-mixing works in 
the EMP stars at AGB phase. Carbon and the s-process 
elements are dredged up to their surface by this mecha- 
nism. 

Observationally, some peculiarities have been pointed 
out, too. It is known that the large portion (20 — 25%) of 
EMP stars show carbon enhancement (Rossi et al. 1999; 
Christlieb 2003) in contrast with a small fraction {^ 1%) 
of CH stars among Pop. II stars. In addition, some 
carbon-enhanced EMP stars stand without s-process ele- 
ment enhancement (Aoki et al. 2007), while no CH stars 
do. The number of EMP stars is very small compared 
to Pop. I and Pop. II stars, which may be referred to 
as the G-dwarf problem in the Galactic halo. All globu- 
lar clusters have the metallicity [Fe/H] > —2.5 and the 
local dwarf galaxies show the metallicity distribution of 
stars that steeply decreases around [Fe/H] ~ —2.5 with 
decreasing metallicity (Helmi et al. 2006). 

In Paper I, we show that the IMF of the EMP pop- 
ulation is different from Pop. II stars, based on the ob- 
served statistics of the CEMP survivors. Only massive 



Komiya et al. 



IMF with medium m.ass Mmd ^ 10 Mq can form sim.ul- 
tancously a large fraction of CEMP stars and the number 
ratio between the CEMP stars with and without the en- 
richment of the s-process elements as observed among 
the EMP stars. This involves the consequence that only 
a small fraction of EMP population stars can be observed 
today as EMP survivors and that most survivors are the 
secondary components of binaries. The effects of high- 
mass IMF on the chemical enrichment are discussed by 
a simple one-zone model in Paper II and shown to be 
consistent with the observed small number of EMP sur- 
vivors. This means that the formation process of EMP 
stars differs from metal-richer Pop. II stars. 

2.2. Metallicity Gap and Hierarchical Chemical 
Evolution Scenario 

Scarcity of stars with [Fe/H] < —4 and the existence 
of HMP stars have drawn wide attention, but are yet 
to be fully clarified. Prantzos (2003) discusses the low- 
metallicity tail of the MDF, suggesting that the prob- 
lem can be alleviated by introducing an early gas infall. 
Karlsson (2005) points out that such MDF break can 
be interpreted as a result of stochastic metal pollution 
process. Karlsson (2006) proposes that HMP stars are 
formed of the matter, pre-poUuted by SNe ejecta and 
polluted by carbon-rich stellar wind following the sup- 
pression of star formation after the first SNe. Since the 
Pop. Ill stars are thought to be massive, some studies as- 
sume the change of IMF into that which enables the for- 
mation of low-mass stars somewhere below [Fe/H] ~ —4 
(e.g., Tumlinson 2006; Salvadori et al. 2006), but some re- 
cent studies show low-mass stars can be formed in a more 
metal poor environment(Omukai et al. 2005; Schneider 
2006). 

We have proposed a new scenario to explain this break 
from the viewpoint of the structure formation process 
in Paper II. According to the current A cold dark mat- 
ter (CDM) model, galaxies are formed hierarchically; 
Low-mass structures are first formed, and then, merge 
with each other and/or accrete matter to form galax- 
ies like the Milky Way. Stellar halos are thought to 
be formed through the accretion of many small halos 
(Searle & Zinn 1978). In early stages of the Galaxy 
formation, star forming gas is incorporated into many 
small halos. The mass of the first star forming mini- 
halos is thought to be ^ 10^ Mq with a baryonic content 
of ~ 2x10^ Mq (Tegmarket al. 1997; Nishi & Susa 1999). 
When a massive first star explodes as core collapse SN, it 
enriches the host mini-halo with iron ejecta of ~ 0.07 M© 
up to 



[Fe/H] = log 



0.07 
2 X 105 



log Afo,0 ^ -3.5, (1) 



on average, where Ape,© is the solar iron abundance. We 
refer to this event as "first pollution", which sets a lower 
limit to the metallicity of the second and subsequent gen- 
erations of stars. In some mini-halos, the metallicity de- 
creases by the accretion of gas and/or by merging with 
mini-halo with pristine abundance after the first pollu- 
tion, which enables the formation of stars with smaller 

metaUicity of [Fe/H] 4. 

In this scenario, stars of metallicity much lower than 
[Fe/H] ~ —4 such as HAIP stars must be assigned to 
the survivors of stars formed prior to the first pollution. 



They are the "first low-mass stars" and their surface 
metal elements are provided either by the surface pol- 
lution after their formation and/or by pre-poUution of 
IGM before the formation of the host halo. We explore 
this scenario taking into account the merging history of 
the Galaxy, and investigate the current appearance of 
the first low-mass stars in Paper HI. In this paper, we 
investigate some models of the metal abundance of HMP 
stars quantitatively and discuss their formation in detail. 

3. MODELS AND COMPUTATION METHOD 



3.1. Merging History 



the Extended Press- Schechter 
Theory 

First, we build merger trees of our Galaxy using the 
extended Press-Schechter (EPS) approach (Bond et al. 
1991; Lacey & Cole 1993). The Press-Schechter theory 
(Press & Schechter 1974) is an analytical formalism that 
can draw halo mass distribution in the universe. The 
extended Press-Schechter theory gives the distribution 
of mass. Ml, of building blocks at redshift zi that yields 
a halo of a prescribed mass M2 at Z2, and provides a 
method of tracing back through the formation history of 
halos. The mass distribution of such building blocks is 
given in the following equation. 



J{Si,uji,S2,uj2)dSi 



UJl — (xl2 



(27r)i/2(5, -52)3/2 



X exp 



{UJI - UJ2Y 

'2{Si-S2) 



dSi, (2) 



where wm 2} = '-^{^{1.2}) is critical density to collapse 
for linear growth theory of overdensity at Z{i,2} s-nd 
'5'{i.2} = S{Mfi^^i2}) is the variance of the density field 
smoothed on mass scale M/j {1 2}- We use this equation 
to determine the distribution of mass Mh^i of the halos 
at zi that eventually form the halos of mass M^^2 at Z2- 

To calculate the merger tree from the EPS theory, 
we adopt the N-fold accretion method constructed in 
Somerville & Kolatt (1999). In this method, a merger 
tree starting from the halo of a prescribed mass branches 
and accretes material as we step back to higher redshift. 

The low mass limit of mini-halos in which stars are 
formed is a critical parameter to plant the merger tree 
and to determine the start of chemical evolution. When 
a halo and gas cloud collapse, the cloud virializes and is 
heated to the virial temperature of Tvir- Dependence on 
mass and redshift at z ^ 1 is described as 

T,i, (xM^/'(l + z). (3) 

It is shown that for rvir(M/i,z) > lO^K, H2 molecules 
form and function effectively as a coolant and the pri- 
mordial gas clouds can cool in Hubble time (Tegmark et 
al. 1997; Nishi & Susa 1999). This sets the lower limit, 
Mfi.i, to the total (baryon -f dark matter) mass of merger 
tree, which is written in the form; 



Mhi^l.Ux 10^ Mg 



1 + z ^ 



Hn 



10 0.6 



100 km s ^ Mpc"^ 

-3/2 /o A \-l/2 



n- 187r2 



(4) 



Here /i is the mean molecular weight of gas: Hq is the 
Hubble parameter: 17^ = rim/i^m + f^fe/(l + z)"^ + 
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Fig. 1. — An example of merger trees for a halo of mass 10 Mq. 
The variations of total (dark matter plus baryon) mass, Mf^, of mini 
halos are plotted against the redshift, z, (bottom) or the age, t, of 
the Universe (top) for 1000 branches randomly selected among the 
total of ~ 2 X 10^. 

^a/(1 + z)^)'- and Ac is the overdensity of virialized 
objects relative to the critical density. As z increases, 
Mhj decreases, i.e., a halo resolves into mini-halos hav- 
ing smaller masses. The low-mass halos at tips of the 
branches of the merger tree are regarded as the hosts of 
the first stars. 

Figure 1 illustrates the merger tree of a halo with 
mass similar to the Galaxy (10^^ Af©) in ACDM universe 
{flM = 1 - f^A = 0.3, Qb = 0.045, h ^ 0.7, as = 0.9). The 
number of branches amounts to ^ 2 x 10^ but here we 
plot only one thousandth of total branches. The first 
mini-halo is formed at 2; > 30 with mass lower than 
10^ Mq. For a merger tree with cosmological parame- 
ters by WMAP 5 years, the number of branches is ~ 10% 
smaller and their typical age of formation is ~ 20% later, 
but our conclusion is almost the same for this model. 
Most mini-halos are formed between z = 20 and z — 5, 
which corresponds to the ages of t ^ 0.15 Gyr and ^ 1 
Gyr from Big Bang, respectively, while some young halos 
are formed at later ages. 

3.2. Star Formation and Chemical Evolution 

We first summarize the assumptions about star for- 
mation of EMP and Pop. II stars and about chemical 
evolution. As for the formation of first-generation stars 
including Pop. HI stars and their feedback effects, we 
test various models as described in the next subsection. 

• Star formation rate (SFR), -ip, in the halo with gas 



M, 



is assumed to take the form of -0 = 



mass, 

c*-/Wgas after the first SNe in host halos, where c* is 
the star formation efficiency (SFE). The value of c* 
is set at 10~^° yr"^ for our fiducial model, which 
means that gas is converged to stars in 10 Gyr. 
We also adopt that c* = 10~^andlO~^^ yr^^ to see 
the dependence on SFE. For the fiducial model, 
SFR becomes ip = lO"** Mq yr~^ for halos with gas 
mass 10^ M0. We register all the individual stars 
for the EMP population and follow their evolution. 

• The mass of each star is specified randomly accord- 
ing to the statistical weight with the initial mass 
function (IMF) at its birth. We assume a lognor- 
mal IMF. 



^(m) oc — cxp 

771 



{log(m/M^d)}' 
2 X Ai, 



(5) 



For EMP population stars, we assume high-mass 
IMF, derived by Komiya et al. (2007), i.e., lognor- 
mal IMF with the medium mass M^^d = 10 Mq and 
the dispersion Am — 0.4. For Pop. II stars with 
the metallicity of [Fe/H] > —2.5, we calculate two 
cases; the same high-mass IMF as the EMP stars, 
and the low-mass IMF, base on the observations of 
the Galactic halo stars of [Fe/H] > — 1.5 (Chabrier 
2003), the latter of which is given by a lognormal 
mass function in eq. (5) with Mmd = 0.22 M© and 
Am = 0.33 for m < 0.8 Mq and the Salpeter mass 
function for m > 0.8 Mq, respectively. 

• Half of all stellar systems are taken to be binaries. 
The distribution, n{q), of mass ratio, q (= TO2/TO1), 
between the secondary and primary stars of binary 
systems is assumed to obey a flat distribution, i.e., 

niq) = 1. (6) 

• Massive stars with to = 10 — 40 Mq explode as core 
collapse SNe after r = 10^ yrs from their birth and 
eject iron yield of mass Mpo = 0.07 Mq on average. 
Very massive stars with to, — 140 — 260 Mq explode 
as PISN and eject iron yield of mass Mpo — 10 Mq 
in the lifetime of r = 2 x 10^ yrs 

• SNe ejecta is assumed to spread instantaneously 
and homogeneously in their host halos. 

• In the flducial model, gas is regarded as retained 
within the host mini-halo. We also compute models 
with gas ejection from the halos. All gas and metal 
are assumed to be blown off the host mini-halo 
when the kinetic energy of a shell formed by super- 
nova is larger than binding energy Ei,in of the gas 
in host mini-halo, i.e., e^EsN > Eurn where Esn 
is the energy ejected by supernova and ek is the 
fraction of kinetic energy. We take Esn — lO^^ert/ 
for type II SN and Esn = lO^^erg for PISN, and 
ek = 0.1 for model with blowout. For simplicity, 
ejected gas and metal are assumed to be mixed with 
gas throughout the Galaxy instantaneously and ho- 
mogeneously. 

• Baryon fraction in primordial mini-halos is as- 
sumed to be the same as that of the average in 
the universe, r^^/^M- When gas and metals are 
blown off by energetic supernova, Mgas vanishes, 
but dark matter halo remains. Afterward, blown- 
off mini-halos accrete dark matter and gas at the 
mass accretion rate described in the merger tree. 
When gas mass reaches Mhji^b/^M, the star for- 
mation starts again. 

Since we are interested in the global features of Galac- 
tic EMP stars, we deal only with iron as a typical indi- 
cator of metal-enrichment in the universe. The detailed 
abundance pattern may depend on the models and yields 
of SN explosions and the dilution process of the ejecta 
in the halo gas. These effects may be subject to large 
uncertainties. 

For metal pollution of IGM, we make simplifled as- 
sumptions and discuss the overall trend. In this paper, 
we apply a semi-analytic method and discuss the depen- 
dence on the feedback effects from massive stars. Galac- 
tic wind from more massive galaxies also contributes to 
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the IGM metal pollution but we neglect them. As shown 
in Section 6, gas blown off the mini-halos by primordial 
stars is the dominant metal source for HMP stars on 
which we are focused. In actually, IGM metal enrich- 
ment is an inhomogeneous process and demands further 
study including effect of the wind originated from Pop. II 
stars (Tornatore et al. 2007; Trenti & ShuU 2010). 

3.3. Primordial Stars and Feedback 

Stars formed without any feedback effect from previ- 
ous stars are referred to as primordial stars. These stars, 
born prior to EMP stars referred to as primordial stars, 
can be very massive and exert various influences on the 
star formation processes of subsequent generations. In 
Table 1, we summarize models parameters for this feed- 
back. We present the most simplified model without any 
feedback from massive stars as fiducial models. In this 
section, we describe other models with feedback effects 
and model parameters. 

3.3.1. Initial Mass Function 

We assume the same high mass IMF with Mmd = 10 
both for EMP stars and for Pop. Ill stars of the fidu- 
cial model, although typical mass of primordial stars can 
be much higher. Previous numerical simulations of the 
first star formation show that, in the first star form- 
ing clouds, only one massive dense clump with mass 
M > 100 Mq is formed (e.g. Bromm et al. 1999; Abel 
et al. 2002). Typical mass, Mmd,p, of primordial stars is 
one free parameter in our model. We calculate models 
with Mmd,p = 10, 100, 200 Mq. For aU the cases, we as- 
sume a lognormal IMF with Am = 0.4. For the model 
with Mmd,p = 200 Mq, most of primordial stars become 
PISNe. For the model with Mmd,p = 100 Af©, stars with 
typical mass end their lives without iron ejection. Re- 
cent simulations show that binary can be formed for first 
stars(Turk et al. 2009; Stacy et al. 2009). In this study, 
the binary fraction and mass ratio distribution of primor- 
dial stars are assumed to be the same as those of EMP 
stars. 

3.3.2. Negative Feedback by Lyman- Werner photon 

For the Pop. Ill star formation, II2 molecules work as 
coolant in mini-halos with Tvir > lO"' K. But once a mas- 
sive star is formed in the primordial gas clouds, whole 
II2 molecules are dissociated by radiation from the mas- 
sive star to quench subsequent star formation (Oniukai 
& Nishi 1999) until it explodes as SN. In addition, the 
formation of low-mass stars may be promoted by SN ex- 
plosions after the first pollution (Machida et al. 2005). 
Therefore, it is speculated that SFE was lower before 
the first pollution. We calculate models with and with- 
out this local negative feedback. For the model with local 
negative feedback, we assume that no stars are formed 
in the mini-halos during the lifetime of first stars. 

We also calculate models with global effect by Lyman 
Werner photons. Lyman Werner photons can dissociate 
II2 molecules also in IGM and suppress star formation. 
We assume that Lyman Werner background from pri- 
mordial stars becomes effective at zi20. We assume that 
all mini-halos with Tvir > lO'^K can form primordial stars 
for z > 20, but in the models with global negative feed- 
back, mini-halos can form stars at 2; < 20 only when ad- 
ditional coolant is available. For more massive halos, on 



the other hand, LW background docs not affect the star 
formation because of self-shielding. Additionally, for all 
the mini-halos with Tvir > W K, we assume that stars 
are formed since gas is ionized by accretion shock and 
cools efficiently (Uehara & Inutsuka 2000). 

The assumption on the LW photons is specified in col- 
umn "Photodissociation" of Table 1. We calculate three 
cases of models without negative feedback, models with 
feedback only in the host halo, and models with feedback 
on the IGM. In the table, these cases are referred to as 
no, local and global, respectively. 

3.3.3. Positive Feedback and Transition of the IMF 

Some possible positive feedback effects on the star for- 
mation from primordial stars are also discussed. They 
can reduce the low mass limit of mini-halos for the star 
formation and also the typical mass of stars. 

One is photo-ionization. Once gas is ionized, II2 
molecules form efficiently with residual free electrons as 
catalysts. Further, HD molecules act as main coolant 
to form stars in the high-density cooled gas. Since the 
first very massive stars emit a large amount of ionization 
photon to ionize the intergalactic medium, some halos 
with Tvir < 10^ K are subject to the effect of the photo- 
ionization before they collapse. In the halos with gas 
once photo-ionized, matter also cools by H2 line emis- 
sions effectively (Ricotti et al. 2002). Ionization can also 
reduce the typical mass of stars of the subsequent gener- 
ations (Uehara & Inutsuka 2000; Yoshida et al. 2008a). 

The other possible catalyst is the metal elements and 
dust. The first very massive stars explode as SNe to pol- 
lute the ambient gas with metals and dust. Line emis- 
sions from metal elements such as carbon and oxygen 
and/or thermal emissions from dust become dominant 
coolants for the polluted gas. Additionally, II2 molecules 
are formed on the surface of dust and lower the typical 
mass of stars formed. Observationally, the detection of 
the HMP stars with [Fe/H] < — 5 indicates that some 
low- mass stars are formed in the gas below [Fe/H] < — 5. 
Bromm & Loeb (2003) propose that carbon of abundance 
larger than [C/H] > —3.5 leads the cloud to fragment 
into smaller clumps. Omukai et al. (2005) argue that for 
gas with metallicity [Fe/H] > —5, the effective cooling 
due to H2 formed on the dust surface enables the low- 
mass fragmentation while gas with [Fe/H] < —6 evolves 
like the primordial gas. Schneider (2006) show that dust 
synthesized by the first supernovae lowers the Jeans mass 
and enables low mass star formation at [Z/H] > —6. 

We investigate some models with the conditions for 
the formation of second generation stars under positive 
feedback. The assumptions are given in column "cool- 
ing" in Table 1. The first model "Tvir > lO^^K" shows no 
positive feedback effect from stars. For this model, stars 
are formed only in halo with Tvir > 10^ K for 2; < 20 
because of negative feedback by LW photons. Models 
with this assumptions are named V. The second model 
"photoionization" shows the effect of photoionization of 
IGM. In the mini-halos formed of prc-ionizcd IGM, gas 
can cool and stars are formed. We treat the ionization of 
ambient gas by the massive stars in the following statis- 
tical manner. One Pop. HI star emits ionizing photons 
and increases the mass, Mion, in the ionized region of 



IGM by AMion = iVion/cscm. 



p, where nip is the mass 



of proton. We use the number of ionization photons. 
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TABLE 1 
Model Parameters 



name 


photodissociation 


cooling 




c*iyr '-) 


MmdA^e) 


<^k 


accretion 


FO 


no 






10-1" 


10 







Fs "^ 


no 






10-10 


10 







F-9 


local'' 






10-9 


10 







F-10 


local 






lO-io 


10 





c 


F-11 


local 






10-11 


10 





c 


FlOO 


local 






10-1" 


100 







FA 


local 






10-11 


10 





A 


FB 


local 






10-11 


10 





B 


FD 


local 






10-11 


10 





D 


Vno 


global^: 


Tvir > 


lO^K 


10-11® 


10® 







V 


global 


Tvi. > 


lO^K 


10-11 


200 


0.1 




Pno 


global 


photoionization 


10-11 


100 







P 


global 


photoionization 


10-11 


200 


0.1 




PIO 


global 


photoionization 


10-11 


10 


0.1 




PlOO 


global 


photoionization 


10-11 


100 


0.1 




M 


global 


[Fc/H] 


> -6 


10-11 


200 


0.1 




M-4 


global 


[Fc/H] 


> -4 


10-11 


200 


0.1 




M-10 


global 


[Fc/H] 


> -6 


10-10 


200 


0.1 





Note. — Model parameters investigated in this paper. Column "photodissociation" denotes the assumption about the suppression of 
star formation by Lyman Werner photon (Sec. 3.3.1). Column "cooling" denotes the criterion for the second star formation with positive 
feedback effects and the change of IMF (Sec. 3.3.2). c* is the star formation efficiency, Mm^.p is the medium mass of the primordial stars, 
and efc is the proportion of SN energy, made up by the kinematic energy for gas blown off the host mini-halos. Last column "accretion" 
denotes the model for the surface pollution of Pop. Ill and EMP stars (Sec. 3.4). See text for detail. 

''For stars with [Fc/H] > -2.5, the low mass IMF (Chabricr 2003) 
is assumed. 

I'ln the mini-halos, star formation is suppressed between the for- 
mation and SN explosion of the first massive stars in these halos. 

'^ Lyman Werner background suppresses star formation of halos 
with Tvir < 10-"' K at z < 20 

TVion, emitted by massive stars by Schaerer (2002) and 
the escape fraction, /esc, from host halo by Yoshida et 
al. (2008a). New halos are formed randomly in the ion- 
ized or unionized region. The probability that a new 
halo is prc-ionizcd in advance of formation is given by 
A^ion/(10"'^^ MQ^^b/^M), in which stars are formed under 
the same condition as the fiducial model. The third is 
the metal-pollution model. In the halos, pre-polluted by 
metals ejected from the mini-halos, gas can cool and form 
stars. We introduce the critical metallicity, [Fe/H]cri for 
the formation of the second generation stars, and test 
the models with [Fe/H]cri = —6 and [Fe/H]cri — —4. 

As stated above, it is thought that primordial stars 
are very massive and EMP stars are less massive than 
primordial ones. There should be a switchover of IMF 
from the very massive (Mmd = 200 M©) to the high- 
mass one (Mmd = 10 Mq) with the formation of low- 
mass secondary survivors in the binaries. We assume 
that all stars formed with positive feedback obey the 
IMF with Mmd = 10 Mq. The IMF of stars, formed 
in the halo with Tvir > lO^K, is also a high-mass one 
with Mmd = 10 Mq. 



white dwarfs with 0.6 M© after their lifetimes of 10"'^'' x 
(M/Mq)-^ yrs. In our assumption, massive stars of 
M > Mup(= 8Mq) explode as SNe, and after that, the 
secondary stars are released from the binary systems and 
become single stars. 

Accreted matter is mixed in the surface convectivc zone 
with mass Mscz- The surface abundance of polluted 
Pop. HI star then becomes 



[Fe/H]p = log 



^Fc,ISM(t)M(t)dt/(MscZ ^Fe,©) 



to 



3.4. Surface Pollution by ISM Accretion 
We assume the Bondi-Hoyle accretion rate, 

M = 47r(GM)V/(V;^ + clf^^, 



(7) 



where p and c^ are the density and sound velocity of the 
ISM, Vr is the velocity of the stars (or the barycenter of 
binaries) relative to the ambient gas, and M is the stellar 
mass (or the total mass of binary). The intermediate- 
mass stars with initial mass M are assumed to become 



(8) 
where to is the time when the star is born, XFo.isM(i) 
is the iron abundance of ISM in their host halo as a 
function of t, taken from our fiducial model with the 
radiative feedback and model with lower SFE (F-IO and 
F-11). Since most observed objects have the mass M ~ 
0.8 Mq, we calculate the surface metallicity distribution 
of Pop. Ill and EMP survivors with 0.8 Mq by assuming 
the surface convection of mass Mscz = 0.003 Mq and 
0.2 Mq for dwarfs and giants, respectively (Fujimoto et 
al. 1995). 

There are several factors that may affect the ISM accre- 
tion and should be considered in the hierarchical struc- 
ture formation. Since the accretion rate in eq.(7) is 
dependent on the dynamics of gas and stellar systems, 
therefore, it may change as halos merge and grow, and 
yet, their effects are not yet revealed. We calculate the 
ISM accretion for four models with simplified assump- 
tions on the dynamics and distribution of stars and gas 
in the mini-halos in the early universe. We use models 
without gas blown off and pre-polluted in order to see 
only the effect of surface pollution to the most metal- 
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poor stars. 

In the simplest model of Case A, we assume that stars 
always move with the virial velocity, i.e., Vvir, in ISM of 
average density equal to the gas density, pvir, of virialized 
mini-halos. Since the virial velocity depends on the halo 

mass [Vvh ex M^' (1 + 2;)^/^] and the density of virialized 
halos is ~ 200 times the mean density of the universe at 
the collapse time [pvir oc (1 + z)^], the accretion rate is 
higher in smaller halos collapsed in higher redshift, as 
given by 



M = 4^(GM)Vvirnb 
-4.3 X 10" 



M 



0.8 Mq 



106 Mp 



1 



10 



3/2 



Mo/yr. (9) 



In actuality, gas and stars are concentrated around the 
center in the mother clouds in which Pop. Ill and EMP 
stars are formed so that the relative velocity is lower 
than Vvir and the gas density is higher than pvir- In 
the primordial clouds, gas cools to ~ 200 K with H2 
molecule as the main coolant and falls toward the center 
to form stars. For the model of Case B, we take Vr to be 
comparable with the sound velocity at this temperature. 
For the density of ISM, we assume an isobaric contraction 
from the virial temperature to 200 K, 



p = pvir (200K/Tv 



(10) 



As the density increases in inverse proportion to the tem- 
perature and the relative velocity decrease as V"^ oc T so 
that the accretion rate increases as M oc pV~^ oc T^^-^. 
Since the masses of host halos of Pop. Ill stars are 
equal or slightly larger than the low-mass limit for star 
formation, the temperature of virialized gas is T ~ 
Tvir{Mfi.i) — 1000 K. Accordingly, the accretion rate may 
be enhanced by 



M oc (200K/1000K) 



-2.5 



56. 



(11) 



After the host halo merges to a larger halo, the stars 
are assumed to move with the virial velocity, K/ir, of the 
merged halo in the ISM of its virial gas density, pvir- 

In addition, we have to take into account the time de- 
lay between the merger of dark-matter halos and the 
change of motion of baryonic systems since the incor- 
poration of over-density region of dark matter as given 
by merger trees does not imply the immediate viraliza- 
tion of the constituent systems. The baryonic systems 
in the merged mini-halo orbit as satellites for some time, 
and the relative velocity of stars to ambient gas may not 
change immediately. Lacey & Cole (1993) calculate the 
time for a satellite's orbit to decay by dynamical friction 
on the basis of the standard Chandrasekhar formula. In 
the model of Case C, we consider this delay time, and 
assume that the relative velocity of stars to ISM remains 
the same as in the satellite baryonic systems for dynam- 
ical friction timescale. Other parameters are taken to be 
the same as in the model of Case B. 

Finally, since most EMP stars are formed as the sec- 
ondary members of binary systems under our high-mass 
IMF, we have to consider the allotment of accreted mat- 
ter onto the binary systems to each member in order to 



estimate the pollution of surface abundances. It is true, 
however, that the gas dynamics of accreted flow in the 
binary systems is yet to be established (Bate & Bonnel 
1997; Ochi et al. 2005). While we simply assume that 
a half of accreted ISM settles onto each member in the 
models of Cases A-C, we compute an additional model 
of Case D where the members of binary systems accrete 
mass independently at a rate given by eq.(7) by taking 
the other parameters in the same way as in the model of 
Case C. 

3.5. Comparisons with Observations 

In our study, we consider not only the pattern of MDF 
but also the total number of EMP survivors observed 
in our Galaxy because of the large survey area and the 
significant efficiency of the identification of stars by the 
HES survey. The effective survey area is 5 = 6726 deg 
(Christlieb et al. 2008). Roughly 40% of the candidates, 
selected by the objective-prism survey, have been ex- 
amined by the spectroscopic follow-up observations with 
medium resolution (Beers et al. 2005). In addition, the 
large limiting magnitude {B < 17.5) assures that almost 
all of the EMP giant stars present in the survey fields 
are identified. Therefore it is expected that the ratio of 
the number of giants detected in the survey to the total 
number of giants in the Galaxy should be S/Aw. As for 
dwarf stars, only nearby stars can be observed because of 
lower luminosity, for which we adopt the observed num- 
ber ratio of EMP dwarfs to giants {— 0.93). 

For stars with [Fe/H] < —2, we use MDF of stars ob- 
served in HES survey for comparison (Beers et al. 2005). 
As for the stars of metallicity below [Fe/H] = —3, it 
is possible only through the follow-up observations with 
high dispersion spectroscopy (HDS) to distinguish their 
metallicity. We use data compiled with SAGA database 
(Suda et al. 2008), which collects all the abundance data 
of metal-poor stars with [Fe/H] < —2.5 obtained with 
the HDS follow-up and published up to date. The follow- 
up HDS observations are reported only for 153 stars for 
[Fe/H] < —3, detected in HES survey and HK survey, 
while HES surveys counts ^ 200 candidates of stars with 
[Fe/H] < -3. For [Fe/H] < -3, we scale MDF of SAGA 
by factor 200/153, since we may well regard the selection 
of target stars for the follow-up observations as being un- 
biased below the metallicity of [Fe/H] ~ —3. In actuality, 
three HMP stars currently known were all identified by 
the follow-up observations of HES samples. We may ex- 
pect to have 3 x (200/153) ~ 4 HMP stars among the 
stars in the fields observed by HES surveys. 

For stars with [Fe/H] > —2, we use kinematically se- 
lected sample of halo stars by Ryan & Norris (1991) 
because the sample of HES survey is biased to low- 
metallicity stars. We scaled their MDF to match the 
MDF by HES survey at [Fe/H] = -2. Observed MDF is 
plotted in Figure 2 by solid line. Recently, Schoerck et 
al. (2008) assess the selection bias of HES survey and 
estimate the corrected MDF. Their result shows a steep 
increase in the number of stars above [Fe/H] > —2.5 as 
plotted in Fig. 2 by dashed line, which is inconsistent 
with MDF by Ryan & Norris (1991), waiting for further 
investigations. 

4. RESULT 1: MDF AND FORMATION HISTORY OF EMP 

STARS 
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Fig. 2. — Basic features of metallicity distribution function 
(MDF) of EMP survivors. Shaded fiistogram denotes the number 
of stars predicted by our fiducial model (FO) with a high-mass IMF 
of the medium mass Mmd = 10 Mq and the dispersion Am = 0.4 
for the same area coverage and efficiency as HES survey. Light 
gray histogram denotes the predicted number when we change the 
IMF into a low-mass IMF of Galactic spheroid (Chabrier 2003) for 
[Fe/H] > -2.5 (Fs). SoUd line shows the observed MDFs for EMP 
stars by the HES survey (Beers et al. 2005) and for the kinomat- 
ically selected Galactic halo stars (Ryan & Norris 1991), matched 
to the HES survey for [Fe/H] > -2.5. For [Fe/H] < -3, plotted 
is MDF of stars with HDS follow-up observation, compiled by the 
SAGA database (Suda et al. 2008). Dashed line shows the "cor- 
rected" MDF of the HES survey by Schoerck et al. (2008). 

In this section, we show results of our fiducial model 
and the models labeled with F without global feedback 
effects. 

4.1. Metallicity distribution function 

We first discuss Model FO; our fiducial model with the 
same high-mass IIVIF and SFE applied to the stars of all 
the metallicity from Pop HI and ElVIP stars to Pop II 
stars of [Fe/H] > —2.5 without any feedback. The resul- 
tant IMDF is shown (shaded histogram) and compared 
with observations in Figure 2. Here the model IVIDF is 
given under the assumption of the same area coverage 
and efficiency as HES survey. In the hierarchical struc- 
ture formation scenario, the stars born before or after the 
first pollution constitute two separate groups. The stars, 
born after the first pollution, reproduce the observed 
MDF in the metalhcity range -4 < [Fe/H] < -2.5. In 
particular, the predicted number of stars agrees well with 
the observed number of EMP survivors. 

Our model well represents the observations in the to- 
tal number of EMP survivors as well as in the shape of 
the MDF. This is due to high mass IMF for EMP pop- 
ulation. The model with low-mass IMF predicts a much 
larger frequency of survivors as shown in the light gray 
histogram in Fig. 2, and is inconsistent with the obser- 
vations. 

We present the results of computations for one merger 
tree because results are almost the same for any merger 
tree. This is because many low-mass mini-halos take 
part in the formation of EMP stars. Merger trees with 
10^2 Mq, have ~ 2 x 10^ branches with 10^ Mq and 
^ 4 X 10^ branches with 10* Mq. The history of a major 
merger in low rcdshift is different from tree to tree but 
the distribution of low mass mini-halos is similar because 
of the large number. For computation with WMAP5 
parameter, merger trees have smaller branches and pre- 
dicted number of Pop. HI stars decreases with ~ 10%. 

In our model, the break of MDF near [Fe/H] ~ —4 
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Fig. 3. — Frequency distribution of metallicity of EMP survivors, 
descendant of SNe in each mini-halo. Thick solid (red) line shows 
MDF of stars, born between the first pollution and the second 
SNe in their host mini-halo. These stars are formed out of gas 
mixed with the ejecta of only one SN and are distributed around 
metallicity slightly larger than [Fe/H] ~ —4. Dashed, dotted and 
thin solid (green, blue, magenta) lines denote the number of EMP 
survivors, descendant of the 2nd, 3rd and 10th SNe, respectively. 

is reproduced without recourse to any ad-hoc change of 
IMF. This is the natural consequence of hierarchical sce- 
nario as discussed in § 2.2. Figure 3 shows the number 
of EMP survivor descendants measured by the number 
of SNe until they acquire their metallicity. The distribu- 
tion of EMP stars, formed after the first SNe and before 
the next SNe in each mini-halos, is plotted in thick solid 
(red) line. The 2nd-generation of stars are distributed 
around [Fe/H] r^ —3.6 and form a steep slope of MDF 
around [Fe/H] ^ —4. Metallicity distributions of sub- 
sequent generations of stars, which are descendants of 
the 2nd, 3rd and 10th SNe, are also shown. The metal 
abundance of EMP stars is approximately proportional 
to the number of progenitor SNe. Our model shows that 

EMP stars with metahicity [Fe/H] 2.5 have ~ 10 SN 

progenitors. 

For [Fe/H] > —2.5, the fiducial model has a similar 
slope of MDF to that of EMP stars since we assume the 
same IMF as for EMP population. This seems consistent 
with the observations for [Fe/H] < —2. Ryan & Norris 
(1991) report that the MDF smoothly extends from the 
peak near [Fe/H] ~ —1.8 into EMP stars for the kine- 
matically selected halo stars. In contrast, the switchover 
to a low-mass IMF, characteristic to Pop. II stars, causes 
a great increase of low-mass survivors to be incompatible 
with the observations, as seen from the Fig. 2. The "cor- 
rected" sample of HES survey (Schoerck et al. 2008) 
has a steeper increase of MDF for [Fe/H] > —3 than the 
sample of Ryan & Norris (1991), and yet, their MDF is 
also by far smaller than the result with low-mass IMF 
and still closer to the fiducial model. This indicates that 
no significant change in the IMF is needed within this 
model as far as the field stars in the Galactic halo are 
concerned. 

4.2. Pop. Ill Stars 

The fiducial model predicts a large number (~ 300) 
of Pop. HI stars. Since there are no other observational 
counterparts in this metallicity range, it is natural to 
assign them to HMP stars. However only three HMP 
stars have been discovered to date and by far fewer even 
if we take into account the fraction of stars with follow-up 
observations. 
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Fig. 4. — MDF resulting from the models under the different 
assumptions on the star formation for Pop. Ill stars. Solid and 
dashed lines denote results of Model F-10 with the negative feed- 
back by the radiation of the first stars and of Model with IMF of 
JW'md = 100 Mq for the stars with Z=0, respectively. 

In Figure 4, we show the model including local radi- 
ation fccdback by solid line. Here we assume that the 
first massive stars of M > 10 Mq, once born, prohibit 
star formation in their host halos until their explosion. 
We find that the predicted number of Pop. Ill survivors 
decreases by a factor of ~ 5. Accordingly, we may expect 
~ 60 Pop. Ill survivors in the HES survey area, which is 
still larger up to by a factor of ^15 than the number of 
currently estimated HMP stars. 

A hump around [Fe/H] ~ —3.5 is due to the star for- 
mation after the first pollution and before the second 
SNe in the host halos, which may be an artifact of our 
assumption to take into account the radiative feedback 
only for the first stars. In this model, the first massive 
star is formed similar to Model FO, but the formation 
of 2nd generation stars, and hence, the second SN ex- 
plosion is delayed. During the elongated time interval 
between the first and second SNe, the second generation 
stars are formed as many as ^ 5 times larger compared 
with Model FO. In actuality, however, radiative feedback 
should work for the second and subsequent generations 
of stars, and hence, the hump will spread to be alleviated 
if we include this effect. 

For models with smaller cg by WMAP 5-year, the pre- 
dicted number of Pop. HI stars becomes ~ 10% smaller 
because of a smaller number of merger tree branches, but 
it is still much larger than the observed number. 

The MDF predicted by the Model FlOO with very mas- 
sive Pop. HI stars is also illustrated in this figure (dot- 
ted line). For Pop. HI stars, we assume the IMF with 
Mynd = 100 M0, and the same set of parameters as the 
fiducial model without radiation feedback. Since we as- 
sume n{q) — 1, the fraction of low- mass star companions 
decreases in inverse proportion to the medium mass. In 
addition, the total number of Pop. HI stars is also re- 
duced because SFE is defined as the total mass of stars 
formed in unit time in unit mass. As a result, the pre- 
dicted number of Pop. HI survivors is reduced by factor 
of ~ 1/100 to be comparable to the number of observed 
HMP stars. However, in this case, significant fraction 
of stars becomes PISNe and the assumption without gas 
blown off the mini-halos is not realistic. We will discuss 
feedback effects with gas blown-off the mini-halos in § 6. 
In this case, the number of stars with —4 < [Fe/H] < —3 
becomes larger than the fiducial model, too, because of 




[Fe/H] 

Fig. 5. — The dependence of the MDF on star formation effi- 
ciency (SFE) for the model with radiative feedback. Dashed, solid, 
and dotted lines show the results for Model F-9, F-10 and F-11 
with c* = 10~^, 10"^" and lO^^^yr^^, respectively. There is a re- 
markable hump around [Fc/H] ~ —3.5 in the model with high SFE 
because the number of the second generation stars formed between 
the first pollution and the second SNe in their host halo increases 
with the star formation rate (oc c*). Difference in MDF reduces 
for higher metallicity. 

a smaller number of Pop. Ill stars than in the fiducial 
model. It is also possible to explain the scarcity of HMP 
stars in terms of lower binary fraction or different mass 
ratio distribution. The formation of the first low-mass 
stars and the origin of HMP stars will be discussed again 
in the following sections. 

4.3. Dependence on the Star Formation Efficiency 

SFE is an important parameter for star formation and 
chemical evolution and yet to be properly revealed in 
the very high redshift and/or in very metal poor envi- 
ronment. It could be different from SFE in the nearby 
universe and is free parameter in our computations. Fig- 
ure 5 compares the results of Models F-9, F-10 and F-11 
with different SFE, c* = 10"^ 10"^° and 10"", respec- 
tively. For the models of two larger SFE, a bump on 
the MDF develops around [Fc/H] ~ —3.5 in proportion 
to SFE, as stated above. For the model of the lowest 
SFE, F-11, on the other hand, no bump is discernible 
since the interval between the successive star formation 
grows as large as the lifetime of massive stars, r ~ 10^ 
yr. Apart from the bump, the MDFs are almost indepen- 
dent of SFE not only for Pop. Ill stars but also for the 
stars of subsequent generations. The reason for this is as 
follows; if we increase SFE, the number of stars increases 
at every timestep in proportion to SFE. However, at the 
same time, the interval of SNe, and hence, the timescale 
of chemical enrichment become smaller. This will reduce 
the number of stars at given metallicity in inverse pro- 
portion to SFE. 

In these computations, we assume that radiative feed- 
back works only for Pop. Ill stars. After the first pollu- 
tion, however, massive stars also dissociate H2 molecule 
and ionize gas in the host halo to decrease the SFE. If 
we assume that one massive star prohibits star formation 
in H II region with mass Mhii, the star formation dies 
down when the prohibited region covers the whole gas 
mass of the host halo, i.e., 






Mmi > M„ 



(12) 



where tms is the lifetime of H II region, /ms is the num- 
ber fraction of massive stars, and M is the averaged mass 
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Fig. 6. — Chemical evolution along the merger tree for Model F- 
11. Most of the EMP stars are formed in high redshift oi z > 10. 
There is a large diversity of metallicity of halos in the high redshift 
and EMP stars are formed in many mini-halos that collapse in 

of'^stars? xiLus the effective star formation efficiency is 
thought to be regulated to 
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under the high mass IMF with M^d = 10 M©. This is 
slightly smaller than what we assumed in our fiducial 
model and seems to be consistent with observed MDF 
without hump around [Fe/H] ~ —3.5. It is true, on the 
other hand, that some positive feedback from massive 
stars is thought to enhance the star formation. In the 
relic site of H II region, the formation of H2 molecule 
is argued to increase the formation rate (Ricotti et al. 
2002). The possibility is explored that SN blast wave 
shell may trigger the subsequent star formation, par- 
ticularly of low-mass stars (Uehara & Inutsuka 2000; 
Machida et al. 2005). These possible effects of positive 
feedback have to be considered in further discussion of 
star formation efficiency in the very beginning of chem- 
ical evolution. In any case, the scarcity of HMP stars 
and the absence of a hump near [Fe/H] ~ —3.5 in the 
observed MDF indicate smaller star formation efficiency 
than Ci, = 10^ ^"yr"^ assumed in our fiducial model for 
the first and subsequent star formation. 

In the following, we show results of models with d, 
10-^^ 



-11 



4.4. Chemical Evolution and Formation History of 
Galactic Stellar Halo 

Figure 6 illustrates the evolution of metallicity of sam- 
pled mini-halos against the redshift and age of the uni- 
verse for Model F-11. Most mini-halos are formed in 
20 > > 5 and they evolve to [Fe/H] > —2.5 in timescale 
for '-^ 10 SN explosions, as shown in Fig. 3. When we ne- 
glect accretion of intcr-halo matter, the timescale, ionri, 
to enrich the mini-halos with metal up to [Fe/H] is given 

by 

C^Mpe/SN 

-68x lOl^'^/^lGyr, (14) 

for Model F-11, where /sn is the number fraction of 
massive stars to be SNe, and M is the averaged stellar 




Fig. 7. — Distribution of the formation redshift of EMP stars 
with [Fe/H] < —2.5 and the mass of their host halos at their birth 
for Models F-11 (top panel) and F-10 (bottom panel), respectively. 
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Fig. 8. — Evolution of the MDF with the redshift. Black, gray, 
cross-hatched, and hatched histograms denotes the MDFs of EMP 
and Pop. Ill survivors formed before z = 20, 15, 10, and 5 , re- 
spectively. Thin solid line denotes MDF at z=0. The number of 
stars in the MDF increases with time, while its shape is almost 
unchanged for [Fe/H] < —2.5 . 

mass. When this timescale is shorter than the collapse 
timescale, icoii, of the mini-halos that enable the star for- 
mation, the chemical evolution progresses separately in 
each mini-halo. For Model F-10 with higher SFE, the 
timescale is short enough for the independent evolution 
and for Model with F-11, some EMP stars are formed in 
the merged halo. 

Figure 7 shows the distribution of EMP and Pop. HI 
stars on the diagram of the formation redshift and mass 
of host halo at their birth for Models F-11 and F -10. 
Most EMP stars are formed in mini-halos with 10® — 
10® Mq at z -- 5-15. As shown in Fig. 6, for Model F-11, 
metallicity z = are 0.7 dex lower than solar metallicity. 
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indicating that SFR is higher at least in the Galactic disk 
at lower redshift. 

We assume that all halos with Tvir > 10^ host stars but 
the re-ionization is thought to prohibit the accretion of 
gas onto halos with small mass of Tvir < iC- For models 
with higher SEE, the results about EMP stars are quite 
similar because most EMP stars are formed at z > 10 
and before the reionization. Eor Model E-11, a significant 
fraction of EMP stars are formed in lower redshift but 
the resultant MDF is quite similar even under the effect 
of reionization. 

Eigure 8 shows the distributions of Pop. Ill and EMP 
survivors that have been formed by certain period or 
redshift for Model F-11. As the redshift decreases from 
z = 20 to 5, the distribution of stars, formed by the 
given redshift, increases greatly in number and extends 
into higher metallicity. Pop. Ill stars and EMP stars 
are mainly formed at high redshifts of z > 10 and the 
shape of MDF remains little unchanged afterwards. In 
the mini-halos, stars have the metallicity of [Fe/H] > —4 
apart from Pop. Ill stars. In the mini-halos that evolved 
to larger mass by accretion and/or by merger, metals 
ejected by SN are diluted in larger mass, and hence, the 
2nd generation of stars tend to have smaller metallicity. 

Consequently, the metallicity of EMP stars weakly de- 
pends on the formation epoch for the hierarchical sce- 
nario. There is a weak age-metallicity trend for EMP 
stars. The trend appears only when tonri([Fe/H]) > icon- 
As seen in eq. (14), the timescale, tenri, of the chemical 
enrichment becomes shorter for higher SEE, and the age- 
metallicity trend among EMP stars is not discernible for 
models with higher SFE. 

5. RESULT 2: ACCRETION OF INTERSTELLAR MATTER 

In our models, a group of stars are formed with Z = 0. 
In this section, we show the results with Models F-11, 
FA, FB and ED, and discuss the surface pollution of 
these stars by the accretion of ISM. In the previous sec- 
tion, we calculate the formation history of metal-poor 
stars and the chemical enrichment history of mini-halos 
with a merger tree. Here we trace the changes of sur- 
face iron abundances through the accretion of ISM for 
all the individual Pop. Ill and EMP survivors formed in 
the model calculation by taking into account the struc- 
tural and chemical evolution of halos. We investigate the 
dependence on the velocity of stars relative to the ISM, 
the distribution of gas in the mini-halos, and the accre- 
tion process in binary systems. Based on the results, we 
may estimate the observed abundances of Pop. Ill sur- 
vivors and discuss the lowest surface metal abundances 
of halo stars. 

For HMP stars, the most metal-poor stars, known to 
date, the peculiar abundance patterns such as very large 
carbon enhancement have been discussed. In this pa- 
per, we focus on the iron abundances, however, since the 
observed abundance patterns of the light elements from 
lithium, carbon through aluminum are investigated in 
detail and shown to be affected by binary mass transfer 
(Suda et al. 2004; Nishimura et al. 2009; Suda & Fuji- 
moto 2009, see § 7.1 below). The iron abundance is a 
good indicator for ISM accretion because it is not af- 
fected by the binary mass transfer and the evolution of 
low-mass Pop. Ill survivor themselves. 



5.1. Surface Pollution History 

We illustrate the detailed history of ISM accretion and 
surface pollution for the Case C model in Figure 9, where 
the time variations in the accretion rates of ISM gas and 
iron and variations in the surface abundances with the 
chemical enrichment of their host halos for three sample 
Pop. Ill survivors in the binaries. In the birth halos in 
which Pop. Ill stars are formed, the ISM accretion rate is 
as large as M ~ iQ-io- -12 Mg/yr. When the birth halos 
merge to larger ones, the accretion rate decreases imme- 
diately by a factor of '--^ 1/100 or less. This is because we 
assume that gas and stars are concentrated around the 
center in their birth halos while the stars come to move 
with the virial velocity after the first merger, which re- 
duces the gas accretion rate by ^ 1/56, as seen in cq.(ll). 
Additionally, the mass of the host halo increases by a fac- 
tor of 2 or more, which further decreases the accretion 
rate, as seen in eq. (9). 

Figure 10 plots the ISM accretion rate against the red- 
shift of their formation for 200 sample Pop. Ill survivors. 
We see that the ISM accretion rate is high for mini-halos 
and it continues until the host halos undergo a merger 
with larger halos. The accretion rate before the first 
merger is a decrease function of formation time of halos 
because of decrease in the pvir and increase in the mass 
of host halos. For a given redshift, it displays a spread 
by a factor of ~ 10 or more, which is caused mainly by 
the different masses of stars or binary systems, since the 
Bondi accretion rate depends on the mass (see below). 

The iron accretion rate increases from zero to 
j^Q-ie 19 M0 yr~^ after the first pollution in the birth 
halos. As the birth halos grow by merger and the ISM 
accretion rate steps down, the iron accretion rate also de- 
creases despite the chemical enrichment of the ISM with 
metals. It decreases typically down to ^ 10^^^ Mq yr~^ 
today. This rate can be changed by formation of the 
galactic disk but the accretion rate today for halo stars 
with high velocity is still much lower than in mini-halos. 
The accretion of ISM and metals is by far more effective 
before the first merger in the earlier universe than in 
larger halos in later time. For Pop. Ill survivors, there- 
fore, the accretion in their birth halos is the predominant 
source of their surface iron abundance. 

Right panels of Fig. 9 show the cumulative masses of 
accreted ISM and iron and the surface iron abundance 
of giant Pop. Ill survivors. Typically, Pop. Ill stars stay 
in their birth clouds in ~ 10^ yr before the first merger. 
During this period, the metallicity of ISM increases from 
Z = to [Fe/H] ^ —2 or more. Pop. HI survivors ac- 
crete ISM to the amount of 1O~^"~^M0 and iron to the 
amount of 10^^" "1° Mq, and which result in the surface 
metallicity of [Fe/H] ~ —4- — 6. After the first merger, 
the accretion rate becomes much lower and the surface 
metallicity little increases. 

The resultant surface abundance strongly depends on 
the merger history of host halos. Between two stars in 
top and bottom panels in Fig. 9, the final surface iron 
abundance differs by up to about 1 dex, though the bi- 
nary parameters are similar. This large difference origi- 
nates in the dependence of the accretion history on the 
formation epoch and mass of birth halos and also on the 
time span to the first merger. 

The accretion and pollution history also depends on 
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Fig. 9. — Surface metal pollution history of three sample Pop. Ill stars of mass 0.8 Mq in the Model F-11. Left panels show the accretion 
rate of iron (thick solid red lines) with the accretion rate of gas (dashed green lines) and the iron abundance of ISM gas (thin solid blue lines) 
as a function of the redshift (bottom axis) and the age of the star (top axis), respectively. Right panels show the surface iron abundances 
of polluted Pop. Ill giants (thin solid blue lines) with the cumulative masses of accreted iron (thick solid red lines) and gas (dashed green 
lines). Vertical lines mark the time for the formation of the stars (thin solid lines), the merger of the birth halos of Pop. Ill stars with larger 
halos (thick solid lines) and the end of the nuclear burning stage of their primaries (dashed lines), respectively. Top, middle, and bottom 
panels represent the results for Pop. Ill survivors, formed with primaries of mass 5.8,29.7, and 6.0 Mq at the redshift z = 17.7,15.1,8.7 
and the time t = 0.20, 0.30, 0.44Gyr, respectively. 



the mass of erstwhile primary stars. Since the Bondi 
accretion rate is in proportion to the mass of systems 
squared, the enhancement of accretion rate in the binary 
is large for low-mass survivors, as given by 

^ l/M^+OSy (15) 

2 V 0.8 y ' ^ ' 

for the low-mass member of 0.8 Mq stars with the pri- 
mary of mass Ml. In the top and bottom panels, we see 
that the accretion rate decreases about by about an or- 
der of magnitude when their primary stars become white 
dwarfs and reduce their mass (vertical dashed lines) . The 
iron accretion rate also decreases, though it is partly 
compensated by the increase in the iron abundance in 



ISM. Since the enhancement of accretion rate ends ei- 
ther when the primary stars lose their masses or when 
the host halos undergo the first merger, there is the op- 
timal mass of intermediate-mass primary stars that can 
bring about the maximum amount of accreted iron: For 
the primary stars that have smaller masses than the opti- 
mal mass, whose lifetimes are longer than the time span 
to the first merger of '^ 10^ yr, the amount of accreted 
iron increases with the mass of primaries as given above. 
For those primary stars that have larger masses and be- 
come white dwarfs before the first merger, the amount 
of accreted iron tends to be smaller because of stronger 
mass dependence of their lifetime. The difference of pri- 
mary mass considerably accounts for the difference in the 
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Fig. 10. — Variation of ISM accretion rate of ~ 200 sampled 
Pop. Ill stars against the rcdshift, the same as dashed hncs in left 
panels of Fig. 9. 

resultant surface iron abundances of these two stars. 

For SN binaries, the accretion rate starts with very 
high rates, as seen in middle panels of Fig. 9. This 
is also the case for the accretion rates exceeding M — 
10-10 M© yr^i in Fig. 10. It has little to do with the 
accretion of iron, however. Since the primary stars are 
likely to be the first massive stars in the halos, no SN is 
expected to pollute ISM of their birth halos until the ex- 
plosion of their primary stars. Before the SN explosion, 
gas is accreted onto Pop. Ill stars but no metal accretes. 
When the primary stars explode as SNe, the secondary 
stars have to be released from the binaries because of 
sudden decrease in the mass of primary stars, and then, 
begin to accrete metals as single stars. 

5.2. Metallicity of polluted Pop. Ill stars 

In this subsection, we discuss the differences of the 
polluted surface abundance in four models of Cases A- 
D. Figure 11 shows the MDFs with the surface metal- 
enrichment due to the ISM accretion for these models. 

For Case A, the surface metallicity of polluted Pop. Ill 
star is [Fe/Hjpoi ^ —6 for main sequence stars and ^ —8 
for giants. The accretion rate is the smallest among the 
four models since the relative velocity between the stars 
and gas in their birth halos is taken to be larger than 
for the other models. Even in this case, the accretion in 
the earlier universe is dominant because of smaller virial 
velocity of birth halos, collapsed earlier. 

For the models of Cases B-D, the ISM accretion en- 
riches the surface of dwarf Pop. Ill stars with metals up 
to [Fe/H] ~ —3. When evolving to giants, the accreted 

metals are diluted to [Fe/H] ^ —5 6, comparable to 

the metallicity observed for HMP stars. Two of three 
known HMP stars are giants and we may regard them 
as Pop. HI stars with the surface pollution. For a sub- 
giant HE1327-2326, the surface convection is shallow and 
likely to be occupied by the matter transferred from the 
primary in which the accreted iron is diluted in the AGB 
envelope. 

The comparison of the models between Cases A and B 
reveals the dependence of accretion rate on the concen- 
tration of gas and on the velocity of stars in the birth 
halos of Pop. HI stars. As shown in eq.(ll), the surface 
pollution rate differs by log (200K/1000K)-2-5 _ 17 jex. 
Because the accretion in the birth halos is dominant, the 
final metallicity also differs by ~ 1.7 dex between Case 
A and Case B. 



Time delay of merger process increases the accreted 
irons by some amount. For Case B, the metallicity dis- 
tribution has low-metallicity tail stretching significantly 
beyond [Fe/H] ~ —8 while most stars fall in the metallic- 
ity of [Fe/H] > -8 for Case C. In the models of Case B, 
some mini-halos merge to a larger halo in short timescale 
from their formation, and hence. Pop. HI stars formed 
in such halos of short lives accrete little ISM and remain 
with very low surface metallicity. On the other hand, 
in the model of Case C, all the Pop. HI stars stay in 
the same circumstances and continue to accrete ISM at 
the same rate as in the host halos during the dynamical 
friction timescale to raise the lowest metallicity. 

From the comparison between the results for Cases C 
and D, we see that the accretion rate rather weakly de- 
pends on the fraction of ISM that goes to secondary. As 
discussed in the previous subsection, the enhancement of 
iron accretion in the binary is limited only to the pri- 
mary stars of intermediate and smaller masses; there is 
a optimal mass of intermediate-mass primary stars that 
bring about the largest effects, while the effect decreases 
both for larger and smaller masses. Accordingly, the en- 
hancement of iron accretion in the binaries ranges from 
several to several tens, but is not so large on average. 

In summary, the accretion of ISM can pollute Pop. HI 
giant stars to [Fe/H] ~ —5 or more if the concentration of 
gas and stars in the central region of mini-halos is taken 
into account. HMP stars can be polluted Pop. HI sur- 
vivors whose surface abundances of iron group elements 
are strongly affected by ISM accretion. The surface pol- 
lution is dominated over by the accretion in the birth 
mini-halos in the early universe, and this indicates that 
we have to take into account the hierarchical structure 
formation history in order to estimate the effects of sur- 
face pollution. Some dwarf stars with —4 < [Fe/H] < —3 
also can be polluted Pop. HI stars but the number of 
those is much smaller than EMP stars with pristine met- 
als. 

Finally, we show the surface metallicity and formation 
redshift of polluted Pop. HI and EMP stars for the model 
of Case C and for the same model but for c* — lO^^^yr^^ 
instead of c* = lO-^^^yr-^ in Figure 12. The surface pol- 
lution rate is almost independent of the formation red- 
shift, and in proportion to SEE. In addition, there is lit- 
tle difference in the formation redshift between Pop. HI 
stars and EMP stars, in particular for the model of larger 
SEE of = 10-i°yr-i: for the model of lower SEE of 
= lO^iiyr^i, the formation of EMP stars is slightly de- 
layed. 

6. RESULT 3: FEEDBACK EFFECTS AND 
PRE-POLLUTION OF IGM 

Primordial stars, formed prior to EMP stars, are 
thought to be very massive and to have global feedback 
effects on subsequent star formation. Metal abundance 
in the most metal-deficient stars can originate in the pre- 
poUution of IGM by these stars. In this section, we show 
the results of models with feedback from primordial stars 
and discuss the effect of pre-pollution on HMP and EMP 
stars. 

Table 2 summarizes model results for the pre-pollution 
and formation of HMP stars. Observationally, 3 HMP 
stars are discovered and several more HMP stars are ex- 
pected in HES survey field if we consider the fraction 
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Fig. 11. — MDFs of Pop. Ill and EMP stars after the surface pollution due to the ISM accretion (solid lines) in comparison with those 
without the ISM accretion (dashed lines). Dark and light gray histograms denote the distributions of polluted Pop. Ill stars for dwarfs and 
giants, respectively. Four panels compare the models computed under the different assumptions on the dynamics of merger and accretion 
processes. Case A: Stars are assumed to move with the virial velocity in the ISM of the virialized gas density of halos (top-left panel). 
Case B: Stars are assumed to move in the cooled and centrally concentrated gas clouds at the velocity equal to the sound velocity of gas in 
their mother halos (top-right panel). Case C: Delay of merging of stars and gas in dynamical friction timescale is taken into account in 
Case B model (bottom-left panel): Case D; The members of binary systems are assumed to accrete mass of ISM independently in Case 
C model (bottom-right panel). 

of stars with HDS follow-up observations and the heavy 
pollution effects on main sequence stars. The first col- 
umn, Af4.5, denotes the expected number of stars, formed 
with the metallicity below [Fe/H] — —4.5. The second 
column, [Fc/H]/GM denotes the final IGM metallicity 
by pre-poUution of blow off metals by supernova. The 
third column denotes the fraction of iron, originated from 
PISN ejecta, among all the iron ejected to IGM. Obser- 
vational absence of nucleosynthetic signatures of PISN 
ejecta indicates that metal from PISNe should not be 
the dominant source of metal for pre-pollution. 

6.1. Negative Feedback by Dissociation Photon 

Top and bottom panels of Figure 13 show the resul- 
tant MDFs for Models Vno and V, i.e., the models with 
suppression of star formation by LW background without 
positive feedback. Hatched areas in Fig. 13 and follow- 
ing figures in this section show the number of stars that 
formed before the first pollution in their host halo, and 
shaded histograms denote the MDF of stars formed after 
the first pollution. We refer to these stars formed be- 
fore the first pollution as local first stars in this paper. 
All pristine iron in their interior comes from intergalactic 
matter polluted by energetic supernovae which blowout 
their host halo. 

For these models, the minimum halo mass for star for- 
mation becomes larger after z ~ 20 because of nega- 
tive feedback with LW background. For model Vno, the 
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Note. — Results for pre-pollution of the IGM. Expected num- 
ber, A'^4.5, of stars with [Fe/H] < —4.5 in the survey area of HES 
survey; metallicity, [Fe/HjiGMi of intergalactic matter at z=0; and 
mass fraction, FepisN/FejGMi of iron ejected by PISNe among to- 
tal mass of iron ejected from mini-halos by SNe. 



number of Pop. HI survivors is much smaller than for 
Model FO or F-11 because the number of primordial mini- 
halos with star formation decreases with the increase in 
the minimum halo mass for star formation. As shown 
in § 5, these stars are possibly observed as HMP stars. 
In such larger halos, SN ejecta suffers from the dilution 
with large mass and the metallicity of stars born after 
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Fig. 12. — Distributions of the giant Pop. Ill {+) and EMP {■) 
survivors on the diagram of the formation redshift, z, and surface 
metallicity, [Fe/H], after the surface pollution (upper section), and 
their metallicity distribution (lower section). Top and bottom pan- 
els show the result for the models with the star formation efficiency 
c* = 10"^'^ yr~^ and 10~^^ yr~^, respectively, and asterisks de- 
note the iron abundances of three HMP stars. 

the first pollution becomes [Fe/H] ^ —5. Break of MDF 
at [Fe/H] ^ —4 does not appear. The predicted num- 
ber of HMP stars is consistent with the observations but 
the assumption that gas is restricted in mini-halo even if 
PISN occurs is not realistic for the low-mass halos. 

For Model V, the predicted number of HMP stars re- 
sults is too small because of the very massive IMF of pri- 
mordial stars. In our model, metallicity of IGM becomes 
[Fe/H] > —4.5 at very high redshift (z ~ 20) by ejecta of 
PISNe. This result indicates that some positive feedback 
effect that enables low-mass star formation is required 
to form HMP stars. At z > 20, few Pop. II stars are 
formed, and galactic wind originated from Pop. II stars 
seems to be negligible for HMP stars. 

6.2. Metal cooling model 

The one possible catalyst for gas cooling and the for- 
mation of the second generation stars is metals. Top, 
middle and bottom panels of Figure 14 show the resul- 
tant MDF of Models M, M-4, and M-10, respectively, 
with positive feedback by metals. In these models, we 
assume that very massive stars are formed in the halos 
of the metallicity below [Fe/H]cri, while the same high- 
mass IMF as for the EMP population is apply to the 
halos of larger metallicity. We assume [Fe/H]cri = —6 for 
Models M and M-11, and [Fe/H]cri = -4 for Model M-4. 

The MDFs of stars formed after the first pollutions are 
very similar to model F and well resemble the observa- 
tions. For Models M and M-10, stars with pre-poUution 
are distributed between —6 < [Fe/H] < —4. They pre- 




[Fe/H] 

Fig. 13. — Resultant MDFs for Models Vno and V. Hatched areas 
denote local first stars which formed before the first pollution in 
their host halo. Shaded histograms show distribution of EMP stars 
formed after the first pollution. Thick solid lines denote observed 
MDF. 

diet 9 and 20 HMP stars, respectively. This is compara- 
ble to observations if we take into account the fraction of 
stars with HDS follow-up and the heavy effect of surface 
pollution on main-sequence stars as discussed in § 5. For 
these models, most metals for pre-pollution are ejecta 
of Type II SNe. This is also consistent with the obser- 
vations. However, the number of local first stars with 
—5 < [Fe/H] < —4 is larger as compared with the obser- 
vations and the metallicity break disappears as a results 
of these stars. If we consider the inhomogeneity in the 
metallicity distribution of IGM, however, the metallicity 
of these stars also spreads over a wider range, and this 
model can be compatible with the observations. 

In Model M-4, we can see MDF break at [Fe/H] -- -4. 
But it predicts too small a number of HMP stars. For 
this model, metals ejected by PISNe overwhelm metals 
ejected by Type II SNe in IGM. More than half of stars 
with [Fe/H] < —3.5 are formed with the metal abun- 
dances of IGM, although no stars are identified with the 
abundance pattern characteristic of PISNe yields. Al- 
though we cannot reject this model because of lack of 
statistics. Model M seems to be more plausible. 

As shown in § 4.3, there is a hump of MDF at [Fe/H] ^ 
-3.5 for Model M-10 with high SFE. Resuhs about pre- 
pollution are similar to Model M and the metal pollution 
of IGM is almost independent of SFE, also. This is be- 
cause metal in the IGM is mostly provided by massive 
local first stars in mini-halos before the first pollution. 
Only one massive star is formed before the first pollu- 
tion in each halo because of local negative feedback, and 
after gas is blown off the halo, no stars are formed in the 
dark matter mini-halos until they accumulate a signifi- 
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[Fe/H] 
Fig. 14.— Same as Figure 13 but for Models M, M-4, and M-10 

cant amount of gas again. Therefore, the pre-poUution 
rate by primordial stars is weakly dependent on the SFE. 

6.3. Photoionization model 

The top panel of Figure 15 shows the resultant MDFs 
for Model P with the photo-ionization positive feedback 
for the formation of the second generation stars. The pre- 
dicted number of HMP stars is slightly smaller than ob- 
served because of the pre-pollution. For this model, MDF 
does not decrease so steeply at [Fe/H] ~ —4. About 
a half of metals for IGM pre-pollution are provided 
by PISNe. Stars, formed out of pre-poUuted IGM be- 
fore the first pollution, are distributed around [Fe/H] = 

—3.5 5, as shown by hatched histograms in the figure. 

More than half of stars with [Fe/H] < —3.5 are the local 
first stars. The abundances of these stars should show 
the abundance pattern characteristic of PISNe yield but 
it should be blended with ejecta of Type II SNe. 

More observational studies of the MDF and the abun- 
dance patterns of stars with [Fe/H] < —3.5 are necessary 
to test the validity of this model. Further theoretical 




[Fe/H] 
Fig. 15.— Resultant MDFs of Models P, PIO and PIOO. 

studies are also required for detailed metallicity distri- 
bution of the local first stars since our results indicate 
only typical metallicity of these stars by assuming ho- 
mogeneous abundance for intergalactic matter. 

6.4. Typical mass of primordial stars 

The predicted number of HMP stars depends on the 
assumption for IMF of primordial stars. Middle and bot- 
tom panels of Figure 15 show the results of Models PIO 
and PIOO. For the Model PIG with Mmd,p = 10 M©, 
many HMP stars are formed from the pre-poUuted gas. 
Because no stars are formed in the mini-halos of union- 
ized gas under the influence of LW background, the 
predicted number of stars with [Fe/H] < —4 becomes 
smaller than Model F but is still much larger than the 
observation. These results indicate that the fraction of 
low mass stars among primordial stars should be lower 
than for EMP stars. 

Despite stars of 200 M© becoming PISNe and stars 
of 100 Mq becoming black holes without metal ejection, 
the results of P with M^d = 200 M© and PIOO with 
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Mynd — 100 Mq are very similar. In both models, a sig- 
nificant fraction of primordial stars become PISNe and 
a small number of HMP stars are formed. This suggests 
the difficulty to know the specific value of typical mass 
and the mass distribution of primordial stars from the 
observations of EMP and HMP stars. 



7. DISCUSSION 



-HMP STARS & FIRST LOW-MASS 
STARS 



In the previous sections, we have shown that the first 
low-mass stars, formed before the first pollution, can 
have the surface iron abundance of [Fe/H] '--^ — 5 through 
the ISM accretion and/or through the pre-poUution of 
energetic SNe in the very early universe to be observed 
as HMP stars. In this section, we consider other proper- 
ties of HMP stars and discuss the formation scenario of 
these stars. 

7.1. Abundances of the light elements 

As for the HMP stars, one of the central issues is 
the origin of their peculiar abundance patterns (e.g. 
see Tumlinson 2007b), reported by HDS observations 
(Christlieb et al. 2002; Frebel et al. 2005; Norris et al. 
2007; Bessell et al. 2004; Aoki et al. 2006). All the 
three HMP stars show large carbon enhancement with 
[C/Fe] > +2 and some stars show the enhancement of 
nitrogen and oxygen also. The enrichments of Na, Ma 
and Al are also reported for two of them but with differ- 
ent extents. 

One plausible scenario for the abundances of these 
elements is the wind accretion of the envelope matter 
ejected from AGB companion in binary systems (Suda et 
al. 2004). In particular, lithium depletion of unevolved 
HMP star, HE1327-2326, indicates that its surface is cov- 
ered by the matter processed by the evolved stars (Aoki 
et al. 2006; Tumlinson 2007b).' Nishimura et al. (2009) 
investigate the nucleosynthesis in the intermediate-mass 
stars with Z — during TP-AGB phase and show that 
the abundance patterns of HMP stars can be explained 
by accretion from AGB companion, including the dif- 
ferences among three HMP stars. Surface carbon abun- 
dances of the main sequence secondary stars are changed 
to [C/H] > —3 by the wind accretion in the binary with 
separation smaller than < lOOOAU (Komiya et al. 2007) 

Maynet et al. (2006) argue that stellar wind from mas- 
sive Pop. HI star with rapid rotation brings about abun- 
dance pattern similar to that of HMP stars. For these 
massive primaries, however, the fast wind velocity pre- 
vents the low-mass secondary to accrete sufficient wind 
matter. The scenario with peculiar SN yield, proposed 
by Iwamoto et al. (2005), has the advantage of letting us 
discuss the peculiar abundance pattern for the light ele- 
ments to the iron group elements in one lump, but diffi- 
culties have been pointed out about the lithium depletion 
in HE1327-2326 as well as about the formation of stars 
with high carbon abundance with the limited amount of 
carbon ejected (see e.g., Nishimura et al. 2009). 

7.2. Number of Pop. Ill and EMP survivors 

Observationally, three HMP stars are detected among 
the sample stars of HES surveys and several more stars 
are expected from the frequency of follow-up with high 
dispersion observations. On the other hand, we predict 
300 stars should be in survey volume in HES survey in 



the fiducial model and 60 stars in the models with lo- 
cal radiative feedback. As seen in Figure 11, they dis- 
guise their surface with the accreted ISM of abundance 
[Fe/H] > —4 during the main sequence stars and possi- 
bly are lost among the EMP stars. Then, the number 
of stars with surface metallicity [Fe/H] < —4 shrinks by 
about half. For model with cosmological parameter by 
WMAP 5-year, the predicted number of Pop. HI stars de- 
creases by ~ 10%, and yet, the predicted number is larger 
than that expected from the observations. This means 
that the formation rate of low-mass secondary stars in 
the binary systems is lower for the primordial stars. The 
number of Pop. HI survivors predicted by some models 
with Mmd ^ 100 Mq for Pop. HI stars reduces to be 
comparable with observations. In conclusion, the IMF is 
of higher-mass and/or the low-mass star formation is less 
efficient in the binary systems before the first pollution 
in the mini-halos. 

Our semi-analytic model lacks information on the spa- 
tial distribution of stars. If there is a significant spatial 
metallicity gradient for EMP stars in the Galactic halo, 
it makes the comparison with the observed MDF not so 
reliable. The spatial distribution of EMP stars yet re- 
mains unrevealed. Recently, CaroUo et al. (2007) report 
that the metallicity of outer-halo stars is lower than that 
of inner-halo stars, but most of their sample stars have 
metallicity of [Fe/H] > —2.5. From our results, we rather 
expect only weak radial gradient of metallicity at least 
for EMP stars since there is little difference in the for- 
mation epoch between the first low-mass survivors and 
EMP survivors as they are formed in many mini-halos al- 
most independently of the redshift, as shown in § 4.4 and 
Fig. 12. For the stars with [Fc/H] > -2.5 in the Galac- 
tic halo, on the other hand, Ryan & Norris (1991) derive 
halo MDF with a peak around [Fe/H] '-^ — 1.8 for kine- 
matically selected sample stars. In our fiducial model, 
the metallicity of low-mass survivors formed at z = 10 

reaches [Fc/H] = —1.5 2. This indicates that the 

formation of stars in the nearby Galactic halo has been 
suppressed after z ^ 10. Residual gas may go to form 
the Galactic bulge and/or the Galactic disk components 
through the merging processes. 

7.3. Formation of the HMP stars 

Our results give insights into the formation of the first 
low-mass stars. The abundance patterns of HMP stars 
show that all are the secondary members of binaries 
with moderate separation and probably with low- and 
intermediate-mass stars. The HMP stars, known to date, 
all exhibit large carbon enrichments, which is possible 
only for the close binaries of separations '-^ 10 — 100 AU 
(Suda et al. 2004; Komiya et al. 2007). This suggests 
that the low-mass stars can be formed exclusively in the 
close binary before the first pollution. 

From a theoretical perspective, it is known that the 
temperature and Jeans mass should be higher for the pri- 
mordial gas clouds than for the star forming clouds with 
metals, because of the absence of the effective coolants 
such as metals and dust. However, when the gas density 
becomes higher than ~ 10^° cm~'^, three-body reactions 
work and promote the formation of the H2 molecule to 
lower the Jeans mass (Palla et al. 1983). It is argued that 
H2 formation on the dust surface also works to change 
thermal evolution from the primordial cloud only in the 
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phase of high density phase ^ 10^° cm^'^ (Oniukai et al. 
2005). The radial extension of such high-density regions 
is ~ 10^ AU. Before the first pollution, therefore, low- 
mass fragmentation can occur only in the contracting 
gas clouds, which may lead to the formation of low-mass 
survivors in close binary systems. Recently, Machida et 
al. (2008) investigate the evolution of first star forming 
core and showed that fragmentation may occur in a later 
collapsing phase. Clark et al. (2008) also demonstrate 
the formation of "first stellar cluster" by fragmentation 
at high density. More investigations are required to re- 
veal the conditions for the formation of the first low mass 
stars, and yet, our results suggest that low-mass stars 
can be formed only in high density cloud cores as sec- 
ondary members of binary systems with lower efficiency 
than BMP stars. 

8. CONCLUSIONS 

We have studied star formation and chemical evolu- 
tion in the early universe from the perspective of the 
hierarchical structure formation of the Galaxy by con- 
structing the merging trees on the A cold dark matter 
scenario according to the extended Press-Schechter the- 
ory. In particular, we investigate the formation history of 
EMP and HMP stars and the origin of iron in HMP stars. 
We demonstrate that the hierarchical chemical evolution 
can reproduce the characteristics of MDF and the fre- 
quency of EMP stars in the Galactic halo by applying 
the high-mass initial mass function and the contribution 
of binary, derived for EMP stars in our previous work 
(Komiya et al. 2007, 2009a). 

Our main conclusions are summarized as follows. 

1. In the hierarchical structure formation scenario, 
the metal-poor stars are divided into two groups, 
according to whether they are born before or after 
the first pollution, i.e., the first SNe pollution of 
the gas in each mini-halo. The stars of the former 
group are made of gas with the pristine abundances 
of mini-halos, and the stars of the latter group are 
enriched with metals, produced by their own SNe 
of each halo. The latter stars well reproduce the 
observed metallicity distribution of EMP stars of 
[Fe/H] > —4 and the observational counterparts 
of the former stars are found in the HMP stars of 
[Fe/H] < -4. 

2. While star formation efficiency has few effects ex- 
cept for the stars of the lowest metallicity, the ini- 
tial mass function directly infiuences the frequen- 
cies of low-mass stars, survived to date for given 
metalhcity. For EMP stars of [Fe/H] > -4, the 
high-mass IMF of medium mass ~ 10 Mq with the 
contribution of binaries, derived on the basis of the 
statistics of CEMP and EMP stars by Komiya et al. 
(2007, 2009a), gives the total number of EMP stars, 
compatible with the observations by HES surveys. 
In particular, the low-mass stars observed in the 
Galactic halo today are mostly (more than 90%) 
born as the secondary members of binaries. On 
the other hand, the low-mass IMF such as derived 
for the Galactic spheroid component gives rise to 
the overabundance of EMP survivors by more than 
two orders of magnitude. 



3. For stars born before the first pollution, the 
scarcity of HMP stars suggests that the star for- 
mation, especially of low-mass survivors, is less ef- 
ficient as compared with that for the EMP popu- 
lation of —4 < [Fe/H] < —2.5. The overproduction 
of low-mass survivors still persists by a factor of 
^ 10 even if we take into account the suppression 
of the star formation efficiency due to the local ra- 
diation feedback of the first massive stars. This is 
indicative of higher-mass IMF and/or less efficient 
binary formation in the gas of the lowest metallicity 
than for the EMP populations. 

4. As for higher metallicity of [Fe/H] > —2.5, there is 
no indication of significant changes in the IMF, in 
particular into a low-mass IMF as long as the fields 
stars in the Galactic halo are concerned. From 
the comparison with the MDF derived by Ryan 
& Norris (1991) for kinematically selected sam- 
ples, however, the star formation efficiency declines 
for [Fe/H] > —1.8. This occurs near the redshift 
z ^ 10 in our fiducial model and may be related 
to the formation of the Galactic bulge and/or the 
Galactic discs. 

5. For the first low-mass stars, we have evaluated the 
effect of ISM accretion after birth by taking into ac- 
count the dynamical and chemical evolution of the 
Galaxy. The accretion in the mini-halos, in which 
Pop. HI stars are formed, dominates the surface 
pollution, giving a much larger amount of accreted 
gas and metals than in the evolved massive halos. 
If we take into account the cooling of gas and resul- 
tant concentration of gas and stars into the central 
region of their birth halos, the ISM accretion can 
give rise to the surface pollution sufficient to ex- 
plain the iron abundances observed for HMP stars. 

6. We have investigated the feedback effect from pri- 
mordial stars, born prior to EMP stars, and the 
pre-pollution of intergalactic gas by SNe in the 
first collapsed halos. Primordial supernovae pol- 
lute intergalactic matter to [Fe/H] = —3.2 3.8 

in our models and some EMP survivors with such 
metallicity can be formed with pre-pollution be- 
fore the first pollution. Metal ejecta from Type II 
SNe that occurred in the mini-halos with positive 
feedback also pollute intergalactic matter. In the 
models with significant feedback effects, the ejecta 
from PISN with a characteristic abundance pat- 
tern should blend and can be consistent with ob- 
servational absence of abundance pattern peculiar 
to PISNe. Detection of 3 HMP stars indicates that 
some positive feedback enables the low mass star 
formation below [Fe/H] < —4. 

From the binary scenario, we have revealed that there 
are two possible origins of iron for HMP stars, the sur- 
face pollution of accreted IMF in the mini-halos and the 
pre-pollution of intergalactic gas with metals that ener- 
getic SNe have blown off the host halos. The polluted 
Pop. HI stars with ISM accretion and the pre-poUuted 
HMP stars differ only in the presence of pristine iron 
in their interior. In order to distinguish these two pos- 
sibilities, further investigations are necessary to reveal 
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the effects that the pristine metals of such small abun- 
dances as corresponding to the iron abundance observed 
for HMP stars have on the nucleosynthesis and evolution 
of low-mass and intermediate- mass stars (e.g., see Suda 
& Fujimoto 2009). 

In this paper, we use the simplified assumption about 
the star formation and deal only with iron. In particular, 
we assume the instantaneous and homogeneous mixing 
of SN ejecta in the halos and in the IGM. In actuality, 
however, we should take into account the dependence of 
the expansion of shells on the explosion energy of SNe 
and host halo mass, and the finite timescale of mixing 
(Greif et al. 2008; Tornatore et al. 2007). It little affects 
the present results as long as the host halos are not dis- 
rupted by SNe, however, since the SNe shells involve the 
circumstellar matter comparable to the smallest masses 
in the mini-halos that can nurture the star formations, 
except that it may reduce the number of stars with the 
smallest metallicity of EMP stars, which are formed in 



the halos of larger masses that collapse later. From our 
homogeneous mixing model, we cannot compute detailed 
MDF of stars formed before the first pollution and addi- 
tional investigation is required to distinguish the origin 
of positive feedback. Further elaborations are necessary 
for the interactions between the star formation and SNe 
and the binary formation under the metal-deficient con- 
dition. We may expect to gain a better understanding 
of EMP stars and the evolution of the Galaxy from the 
large-scaled, deep surveys of halo stars now in progress 
such as SDSS/SEGUE and LAMOST, as weU as from 
the detailed abundance studies with large telescopes. 
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